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The development of a working quantum computer utilizing electron spin states as qubits is a 
major goal for many theorists and experimentalists. The future applications of quantum information 
technology would also require a reliable method for the transportation of quantum information. A 
promising such method is the propagation of electrons by a surface acoustic wave (SAW). In this 
paper, we simulate the SAW transfer of two interacting electrons through a Y-shaped junction. 
Our results show that the singlet and triplet states can be differentiated in the Y-junction by an 
electrostatic detuning potential, an effect that could be used in for example measuring the state of 
a two- spin qubit. 

PACS numbers: 73.63.Kv,03.67.-a,73.21.La 



The spins of electrons are promising candidate for a qubit 
A framework for using two-electron spin eigenstates 
as qubits was proposed by Levy in 2002 [2[. After that, 
a lot of progress has been made on this direction Q. 
The coherent electrical control and measurement of two- 
spin qubits has been presented experimentally Q . Micro- 
scopic theory [5] can explain experiments and predict 
new phenomena Q • Remarkably long dephasing times in 
two electron spin states have been shown [8]. Recently, 
entanglement between a pair of two-spin qubits has been 
demonstrated experimentally [sj]. This is an important 
achievement, as entanglement is essential in a working 
quantum computer. 

The two-electron spin eigenstates are the singlet state 
and the triplet states. The different symmetry proper- 
ties of the spatial wavefunctions affect the behavior of 
the electrons. Electrons in the singlet behave effectively 
as bosons. The so called exchange force, a geometrical 
consequence of the symmetrization of the spatial wave 
function, causes an attraction between the electrons in 
the singlet state whenever their wave functions overlap. 
For the triplet state the effect is the opposite, the overlap 
of the wave functions results in an effective repulsive force 
between the electrons. This phenomenon is important for 
the quantum information technology and we show that it 
can be used to differentiate the singlet and triplet states 
or for example to measure the spin state of the system. 

Few-electron circuits of quantum information technol- 
ogy, consisting of networks of quantum dots, require a 
method for transportation of electrons between different 
parts of the circuit. A very promising method for trans- 
ferring electrons between quantum dots is by a surface 
acoustic wave (SAW). A framework for SAW transport 
based quantum computing was presented by Barnes et 
al. [10]. It has been shown experimentally that single 
electrons can be captured from one quantum dot and 
transported into another by a SAW pulse moving along 
a one dimensional channel |l2j. The SAW induced 



In this letter, we simulate a system of two electrons 
moving along a SAW pulse in a Y-shaped channel in 
which the Y-branches are detuned using electrostatic gat- 
ing. The potential minimum created by a SAW pulse can 
be thought as a moving quantum dot that confines the 
electrons during their transfer To|, [Hj . As the SAW pulse 
moves at a constant velocity, the system can be modeled 
as a stationary quantum dot that splits into two separate 
dots at the Y-junction, see Fig. [TJ We show that, due 
to the exchange force, the Y-junction and the detuning 
potential can be used to split the electrons in the triplet 
states into the two different Y-branches, while the elec- 
trons in the singlet state are both transferred into the 
Y-branch that is energetically more preferable due to the 
detuning. 
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transfer of electrons has also been realized in the cases of 
Y- and X-shaped junctions 
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FIG. 1. (Color online) SAW induced transfer of electrons 
through an Y-shaped junction. The potential minima caused 
by the SAW pulse are modeled as elliptic quantum dots mov- 
ing at a constant velocity. The initial dot is split into two 
similar dots when it meets the Y-junction at the time t — t\. 
The singlet state (denoted by S) is transferred to the Y-branch 
that is energetically preferable due to a detuning potential, 
while the triplet state (denoted by two spin-up arrows) is 
split into the two Y-branches. 
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We model the SAW transfer of a two-electron system 
through a Y-j unction with the time-dependent Hamilto- 



m = E 



i=i 



V int (r u r 2 ). (1) 



Here, Vi nt is the pairwise Coulomb interaction, and m* 
is the effective mass of an electron. The magnetic field is 
omitted in our computations. 

The external confinement is initially an elliptical 
quadratic well centered at the origin. At time t = ti, 
the quantum dot is split into two dots that start to move 
apart from each other. The external confinement poten- 
tial V ext (rj, t) is written as 

Vw 2 (amm{\ Xj - a(t)/2\ 2 ,\ Xj + a(t)/2\ 2 } + f3y 2 )+Dxj 

(2) 

Here, a(t) is the distance of the dot minima. a(t) = 
when t < t\ and when t > t\ it is an increasing function 
of t. A linear electrostatic detuning potential of strength 
D is included. With D = 0, the minima of the dots are 
located at points (x,y) = (±a(t)/2, 0). A non-zero D 
shifts the locations of the minima slightly (their distance 
is still a(t)). In Eq. ((2]), a and /3 define the x- and y- 
direction confinements of the dot. A schematic of the 
Y-junction can be seen in Fig. [T] 

The ground state of (pQ) is computed by using the 
exact diagonalization (ED) method. The Fock-Darwin 
(FD) states are used as the one-particle basis in the ED. 
The matrix elements for the Coulomb interaction and the 
parabolic external confinement can be computed analyt- 
ically for the FD-states. The Hamiltonian is stored as 
a sparse matrix and its ground state energy and eigen- 
vector are computed by the Lanczos-iteration. The time 
evolution of the system is computed by propagating the 
initial ground state ^(0). 



tj){t + At) = exp 



--H{t)AtU(t), (3) 



where the matrix exponential is evaluated using the 
Lanczos algorithm. The charge density was computed 
from the two-body wave function using the reduced one- 
particle density matrix. 

The two-body wave function, a product of the spin- 
and spatial parts, is anti-symmetric. A symmetric spin 
part corresponds to an anti-symmetric spatial wave func- 
tion, and vice versa. The four lowest eigenstates of a 
two-electron system are the spatially symmetric singlet- 
state \S), and the spatially anti-symmetric triplet-states 
|T_), | To), and |X+). In the zero magnetic field, the three 
triplet states are degenerate and the singlet is always the 
ground state. In the Lanczos-method, the triplet-states 
can be obtained by fixing the z-component of spin to ±1. 

In our simulations, the confinement strength was set to 
hujQ = 3 meV. Several different shapes of elliptical dots 



were used (the values of a and f3 in Eq. (j2j)). The typical 
GaAs parameter values for the effective mass and the di- 
electric constant are m* = 0.067m e and e r = 12.7. After 
the split at t = ti, a(t) was increased linearly from zero 
to the value of 100 nm. The linear detuning strength 
D was set to be in the same order as the linear part of 
the parabolic confinement, D ~ Do = x 30 nm. 

In experiments, this would correspond to two gates with 
a potential difference of approximately 50 mV set 200 
nm apart from each other. The many-body Hamiltonian 
was created using the first 120 FD-states in energy. The 
convergence of the singlet- and triplet-energies was stud- 
ied with respect to the basis size, and 120 FD-states was 
found to be enough to obtain accurate results up to the 
distance 100 nm between the dot minima. 

We did the simulations using several different time win- 
dows and detuning strengths. The time step At was set 
to At < 1 ps. This was found to be small enough to 
ensure accurate computation of the dynamics. It should 
be noted that the simulation can be started at t = ti, as 
the system is in the ground state at t < t\. The relevant 
timescale in the SAW transfer through our simulated sys- 
tem is from 10 ps to the nanosecond scale [ll|, [l2, 14]. 

The results of our simulation in which a is increased 
form to 100 nm in the time of 1 ns can be seen in Figs. 
[2] and 03 The figures were obtained with approximately 
spherical dots, a = 0.9999 and /3 = 1. We included a 
small radial asymmetry to align the electrons along the 
x-axis in the initial dot. The detune strength is set to 
be D = 0.20 x m*^ x 30nm. The number of charges 
in the left and right dots are shown as a function of a 
in Fig. O The number of electrons in the left dot was 
obtained by integrating numerically the charge density 
in the half plane x < #o, where xo = —D/m*^ is the 
middle point between the dots. As for the right dot, the 
integration was done in the area x > xq. Fig. [3] shows the 
evolution of the electron density during the simulation. 
The densities for both the singlet- and triplet-states are 
shown at 20 nm intervals, starting from a = nm. 

The singlet electron density in Fig. [3] shows initially 
just one large central peak. The triplet density shows two 
very distinct peaks corresponding to the two electrons 
even at a = when there is just one potential minimum. 
This difference is explained by the symmetry properties 
of the many-body wave functions; the anti-symmetry of 
the spatial wave function of the triplet state results in 
the density vanishing when the electrons are near to each 
other, i*i ~ T2- 

As the dots are pulled apart from each other, the 
singlet-electrons start to localize more into the left dot, 
as can be seen in Fig. [2j The density in Fig. [3] shows 
that at first a second peak starts to form in the right 
dot. However, the amplitude of this second peak de- 
creases (compared to the left peak) as the distance of the 
dots increases. The right peak disappears when the dis- 
tance of the dots is about 80 nm. When a = 100 nm, 



3 




20 40 60 80 100 
a (nm) 

FIG. 2. (Color online) Number of electrons in the left- (solid 
lines) and right-dot (dashed lines) as a function of the dot 
distance a. Singlet (blue/dark grey) and triplet (orange/light 
grey) results are both shown in the figure. The detuning 
strength was set to D = 0.20 x m*cJo x 30nm, and the time 
window is 1 ns. 



both electrons are fully localized in the left dot, as can 
be seen in Figs. [2] and [3] 

The triplet state behaves in a qualitatively different 
way as the dots are pulled apart. The density peak in 
the right dot never vanishes, it just moves further right. 
At a = 100 nm, the density goes to zero between the 
dots, and there is one electron in the left dot and one in 
the right. The electron density is symmetric with respect 
to the two dots, as can be seen in the integrated electron 
density in Fig. [2] The spatial antisymmetry and the 
resulting repulsive exchange force does not allow both the 
triplet electrons inhabiting the energetically preferable 
left dot. 

The rate at which a is increased was found to affect 
the results. If a grows slowly, the time evolution is ap- 
proximately adiabatic; at the time t the system is in the 
ground state of Eq. (pQ). The results shown in Figs. [2] 
and [3] (the time window is 1 ns) are very close to the adi- 
abatic behavior of the system. When the time window 
is decreased to the value of 10 ps, the singlet state is no 
longer transferred as completely into the left branch as 
in Fig. [3l A small but visible peak stays in the right 
dot even at a = 100 nm. As the window is further de- 
creased, more and more charge is left in the right dot at 
the end of the simulation. It should be noted that in our 
simulations, the time window length corresponds to both 
the SAW pulse speed and the angle of the Y-j unction. 
The higher the angle (relative to the ?/-axis in Fig. [1]) 
is the faster the distance between the dots at the two 
Y-branches increases. 

The behavior of the singlet- and triplet-states was 
found to be quite sensitive to the strength of the detun- 
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FIG. 3. (Color online) The evolution of the charge density as 
a is increased from nm to 100 nm. The density is shown at 
20 nm intervals, starting from a = 0. The blue plots in the 
left show the density of the singlet state, and the orange plots 
in the right the triplet state(s). The detuning strength was 
set to D — 0.20 x m*w§ x 30nm, and the time window is 1 ns. 
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FIG. 4. (Color online) The effect of the detuning strength D 
on the behavior of the system. In the left: D — O.I6D0, with 
D m*cjg x 30 nm. The upper plot shows the number of 
charges in left and right dots as a function of a (see Fig. [2]). 
The lower plot shows the singlet charge density at a = 100 
nm. In the right: D = 0.23Do- The upper plot shows the 
number of charges a function of a (see Fig. [2]). The lower 
plot shows the triplet charge density at a — 100 nm. The 
results are obtained adiabatically. 

ing. If the detuning was too low, the singlet state would 
also split, and if it were too high, both the singlet and the 
triplet would transfer to the left. The region in which the 
triplet would split and the singlet would not was found 
to be quite narrow, from D = 0.17 Dq to D = 0.22Dq, 
with Do = rn*ujQ x 30nm (a was increased from nm to 
100 nm). Cases of too low and too high detunings are 
shown in Fig. HJ The results are obtained adiabatically. 
In the left plots, the detuning is D = 0.1 6Dq and in the 
right plots it is D = 0.23A)- In the D = 0.16A) case, 
the singlet is split at a = 100 nm, although it seems to 
be transferring initially to the left as in Figs. [2] and [3j 
In the D = 0.23Z}o case, the triplet also localizes to the 
left. 

The effect of the detuning can be analyzed qualita- 
tively by looking at the energy levels of the system. As a 
increases, the dots slide up (the right dot) or down (the 
left dot) the linear slope. The linear detuning lowers the 
potential of the left dot by a value Vd{o) that increases 
with a. In the split state, when one electron inhabits 
each dot, and the distance of the dots is large, the elec- 
trons can be thought as non-interacting. The electrons 
stay split until the detuning lowers the potential energy 
of the left dot enough to overcome the Coulomb repul- 
sion between them, i.e. the distance of the dots become 
high enough. As the triplet is higher in energy than the 
singlet, there exists values of a and D so that the triplet 
is split while the singlet is not. 

However, when a is increased enough, the slope will 



eventually win in both the singlet- and the triplet-cases. 
This might not hold anymore for very large displace- 
ments a though, as the middle ridge may become too 
high for the electrons to tunnel through it. In real ex- 
periments, the dot distance (the maximum distance of 
the Y-branches) is naturally limited to some maximum 
value, and the detuning can be set accordingly. 

The dimensions of the dots, the values of a and f3 in 
Eq. (|2|), also affect the behavior of the system. If f3 < a, 
it is preferable for the electrons to align along the 7/-axis 
(see Fig. [T]) in the initial dot at t < t\, which results in 
both the singlet and the triplet transferring to the left. 
Consequently, only values a < j3 were studied. These 
values affect the energy levels of the system and thus 
also the behavior of the electrons in the Y-junction. If 
a < /?, the detuning could be set so that only the triplet 
would be split at the end of the simulation with a — 100 
nm. 

In summary, we predict that a Y-shaped junction with 
an electrostatic detuning potential can be used to differ- 
entiate the two-body singlet- and triplet states in SAW 
induced electron transport. If the detuning strength 
is suitable, the triplet states are split into the two Y- 
branches due to their repulsive exchange force, while the 
singlet state is preserved and is transferred into the ener- 
getically preferable branch. This phenomenon could be 
used in quantum information technology. For example 
in quantum computing, it could allow measurement and 
control of two-electron qubits. 
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